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Abstract 

As the matter produced in a relativistic heavy ion coUision cools through the 
QCD phase transition, the dynamical evolution of the chiral condensate will 
be driven out of thermal equilibrium. As a prelude to analyzing this evolu- 
tion, and in particular as a prelude to learning how rapid the cooling must 
be in order for significant deviations from equilibrium to develop, we present 
a detailed analysis of the time-evolution of an idealized region of disoriented 
chiral condensate. We set up a Langevin field equation which can describe 
the evolution of these (or more realistic) linear sigma model configurations 
in contact with a heat bath representing the presence of other shorter wave- 
length degrees of freedom. We first analyze the model in equilibrium, paying 
particular attention to subtracting ultraviolet divergent classical terms and 
replacing them by their finite quantum counterparts. We use known results 
from lattice gauge theory and chiral perturbation theory to fix nonuniver- 
sal constants. The result is a theory which is ultraviolet cutoff independent 
and that reproduces quantitatively the expected equilibrium behavior of the 
quantum field theory of pions and a fields over a wide range of tempera- 
tures. Finally, we estimate the viscosity r]{T), which controls the dynamical 
timescale in the Langevin equation, by requiring that the timescale for DCC 
decay agrees with previous calculations. The resulting r]{T) is larger than 
that found perturbatively. We also determine the temperature below which 
the classical field Langevin equation ceases to be a good model for the quan- 
tum field dynamics. 

PACS: ll.lO.Wx, 12.38.Mh, ll.30.Rd, 25.75.-q 
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I. INTRODUCTION 



The relativistic heavy-ion colhder (RHIC) is now creating strongly-interacting matter 
at high enough energy density that there is every expectation that temperatures above 
the QCD phase transition are being explored. The question of equilibration is crucial to 
the understanding of these collisions. There has been much progress on one aspect of 
this question recently: we now know that (at least at arbitrarily high collision energies) 
an equilibrated quark-gluon plasma will be created [Q. Furthermore, the first analyses of 
elliptic flow data can be interpreted as indicating early thermalization in RHIC collisions . 
What remains quite unclear, however, is whether even if an equilibrated partonic medium is 
attained early in the collision, thermal equilibrium is maintained as the matter cools through 
the QCD phase transition and becomes a gas of hadrons, which subsequently freeze out and 
are seen in a detector. 

The QCD phase transition is likely a smooth crossover if it is traversed at low baryon 
number chemical potential fiB- At larger /x^, it is thought to be first order. This means 
that at a critical ij,b, there is a second order phase transition in the universality class of the 
three-dimensional Ising model • The nature of the possible nonequilibrium phenomena 
are quite different in these three cases [0]. If the transition is first order, bubbling may yield a 
spatially inhomogeneous final state . If the transition occurs near the second order critical 
point, then nonequilibrium effects tend to obscure the unique fluctuations characteristic of 
the equilibrium critical point: in equilibrium, the correlation length would diverge there; 
nonequilibrium effects smooth out this divergence 0. The farther from equilibrium the 
evolution is, the less likely are distinctive signatures of the critical point. Finally, if the 
transition is a crossover, as is likely at RHIC where /i_B is small, and if this crossover is 
traversed sufficiently quickly that the dynamics can be treated as a "maximally out of 
equilibrium" quench, then long wavelength oscillations of the chiral condensate are excited 
to amplitudes which are greatly in excess of those present in equilibrium |jTO|. The primary 



signal of these disorientations of the chiral condensate (DCC) is a large number of soft pions 
exhibiting large fluctuations in isospin |11"14,1C], which have been looked for and not seen 
in lower energy heavy ion collisions at the CERN SPS []T3[. Many other DCC signatures 



have also been discussed [|^]. The crucial unresolved question in all three cases is: "How 
rapid must the cooling be in order for the evolution to deviate significantly from thermal 
equilibrium?" In this paper, we lay the ground work for answering this question in the case 
of a smooth crossover. 

We construct a Langevin description of the nonequilibrium dynamics of the chiral con- 
densate which we can use over a range of temperatures around the phase transition. Our 
dynamical degrees of freedom are those of the classical linear sigma model. This is valid as 
long as rricr and m^r are much smaller than the temperature. It is therefore valid near Tc 
sufficiently close to the chiral limit. We shall use it with parameters chosen to reproduce the 
physical pion mass at zero temperature, and this will mean that we are using the classical 
theory beyond its realm of quantitative validity. Until quantum field theoretical approaches 
(for example those of Refs. [|T^) mature to the point where they can be used to describe 
spatially varying field configurations, the present classical analysis is of value as a guide. 

Why introduce a Langevin heat bath? That is, why go beyond the purely classical 



analysis of Refs. |10,18]? Allowing the long wavelength degrees of freedom to exchange 
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energy with a heat bath is a crude way of representing the existence of other degrees of 
freedom. Our reason for introducing it is that we wish to investigate the response of the 
long wavelength degrees of freedom to decreases in the temperature T of the heat bath 
that occur with varying speeds. If T is reduced arbitrarily slowly, the system must stay 
arbitrarily close to equilibrium. If T is reduced suddenly, nonequilibrium dynamics results. 
With a Langevin equation, therefore, we shall be able to learn how fast the cooling must be 
in a heavy ion collision if the traversal of the crossover region of the phase transition is to 
be associated with DCC phenomena. In addition to a time varying heat bath temperature, 
we shall also be able to use the freedom to let T vary as a function of position to consider 
initial conditions in which the hot region is finite in extent. 

We leave the answer to the questions just posed to a subsequent paper. Here, we focus 
on setting up the Langevin description, and choosing all the associated parameters. Many 
parameters can be chosen entirely with reference to properties of the system in thermal 
equilibrium. We can fix these parameters by comparison with known results from chiral 
perturbation theory and lattice gauge theory. In so doing we obtain a Langevin equation 
which gives a reasonable description of the equilibrium physics from very low temperatures 
up to temperatures well above Tc. This is possible in equilibrium, but the classical description 
of dynamics must break down at low temperatures, where mo-,7r 3> T. We discover this 
breakdown when we seek to fix those parameters in the Langevin equation which require 
dynamical input, namely the viscosity 1]{T). To this end, we analyze the decay of an idealized 
DCC with time, due to the presence of the heat bath. This problem has been analyzed 



previously by Steele and Koch ]T9[, and we use their results to fix the viscosity in our 
Langevin equation. We find that in the vicinity of Tc, where we expect our treatment to be 
valid, we can successfully reproduce the DCC decay timescale of Ref. [|l^]. In order to do 
so, the T] we introduce is about a factor of five larger than that predicted by perturbative 
calculation [^. At lower temperatures, the whole analysis breaks down: even with the 
dissipation due to the heat bath completely removed by setting 77 to zero, we find that the 
DCC decays much more quickly than it should. The long wavelength DCC is decaying not 
via interaction with the stochastic heat bath but via its interactions with modes which we 
are treating classically. We are therefore able to evaluate the temperature below which our 
classical Langevin field equation should not be used to model the dynamics. 

Our analysis complements that in several previous papers. Greiner and Miiller ||21| 
and Rischke have derived effective Langevin equations by integrating out hard modes 
perturbatively and found that these typically include correlated and multiphcative sources 
of noise and dissipation. Xu and Greiner then explored the effects of such noise and 
dissipation on + 1-dimensional dynamics. This treatment of the noise goes beyond ours, 
since we treat the heat bath as an uncorrelated source of fluctuation. However, whereas 0+1- 
dimensional calculations (like those of Refs. [^,^) treat the dynamics of only one mode 
explicitly and thus assume spatial homogeneity, our 3 + 1-dimensional Langevin equation 
allows an analysis of more realistic, spatially varying, configurations. The 3 + 1-dimensional 
Langevin equation has been introduced previously by Chaudhuri although as he notes in 
Ref. pS) all of his results depend sensitively on his choice of lattice spacing. We introduce the 
correct T-dependent counterterms in order to obtain results at long wavelengths which are 
independent of the lattice-spacing. Randrup has addressed the linear sigma model dynamics 
in a series of papers, using either the classical theory or a Hartree approximation pB[ . 
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These studies have shown the importance of accounting correctly for the multidimensional 
expansion of the hot system, and have provided initial answers to some of the questions we 
wish to investigate. The stochastic linear sigma model we develop below has the relative 
advantage of having dynamics which equilibrate at late times in the absence of expansion 
and in addition allows us to take some account of the effects of shorter wavelength degrees 
of freedom on the longer wavelength modes we treat classically. 

In Section II we describe the model, paying close attention to the counterterms which 
are needed in order to ensure that the long wavelength physics is not divergent as the lattice 
spacing is taken to zero. In Section III, we evaluate the order parameter as a function of 
temperature and fix the most important finite counterterms by enforcing agreement with 
chiral perturbation theory at low temperature and with lattice gauge theory calculations 
of the critical temperature. We also evaluate the temperature-dependent pion and sigma 
masses. In Section IV, we analyze the decay of a DCC By comparing the DCC lifetime 
we obtain to that calculated previously by other methods, we fix the viscosity r] in the 
Langevin equation, which couples the classical fields to the heat bath. We also determine the 
temperature below which a classical field treatment fails to describe the dynamics correctly. 
We close in Section V with a look towards the future. 



II. LANGEVIN EQUATIONS AND COUNTERTERMS 



At low temperatures, the correct effective field theory for QCD is the nonlinear sigma 
model, describing the dynamics of the pseudoscalar pions. At the QCD phase transition, 
however, chiral symmetry is approximately restored and the scalar sigma becomes approx- 
imately degenerate with the pions. These four degrees of freedom are the lightest in QCD 
near its finite temperature phase transition. Indeed, if the up and down quarks were mass- 
less, these four modes would all be massless at Tc, making the transition an 0(4) second 
order phase transition |^|14|. Lattice QCD calculations suggest that with quark masses 
as in nature, these modes have inverse correlation lengths which are comparable to ||28|| . 
This means that although these modes are still the lightest degrees of freedom, the classical 
treatment we employ is at or beyond the limit of its regime of validity. Two asides should 
be noted at this point. First, at the /z^ at which there is a second order critical point in the 
QCD phase transition, the sigma becomes massless while the pions do not. Second, there 
have been recent suggestions that the lightest scalar glueball may play an interesting role at 
Lattice QCD calculations show that this has a correlation length close 

and we therefore neglect it. 



the transition 

[BOl, much shorter than that of the pion pi 



to (3Tc 

Writing the pions and the sigma as an 0(4) vector 
Lagrangian is given by 



(a, tt), the linear sigma model 



(1) 



with summation over repeated indices implied. The context in which this is of precise validity 
is two-flavor QCD in the chiral limit, on the assumption that this theory has a second order 
phase transition, as is consistent with what we know from lattice QCD calculations P2 . 
Then, that transition is in the same 0(4) universality class as the second order phase tran- 
sition found at nonzero temperature in the theory (|l]) with H = [^,0. The explicit 
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symmetry breaking term H, the coupling A, and the vacuum expectation value v can be 
expressed in terms of the phenomenological parameters m„, m.,^, and at tree level via 



2 2 

mt — m 



ml — Sml 



H = Umt, A = :::^-— v'= "\ ^"- /.^ (2) 

2/2 ml - ml 

We take (/^,m^,m<,) = (92.4,135,600) MeV, meaning H = (119. MeV)^ A = 20.0 and 
V = 87.3 MeV. 

The simplest Langevin prescription corresponds to including additive stochastic sources 
^a(x,t) and dissipation terms in the Euler-Lagrange equations of motion for the fields: 

^ - V^0, + X{<Pl- v') 0, - H6^, = -nJ-^ + . (3) 
The stochastic fields are taken to obey the fluctuation-dissipation relations 

{U^)Uy)) = '^Vai>{x)T5\x-y) , (4) 

with {iai.^)) = 0. In this paper, we shall always choose rjab{x) = rj{T)5ab- This dynamics 
reduces to the classical microcanonical relativistic evolution for the pions and sigma in the 
limit of small dissipation 77 — > 0. In the limit of vanishing masses the long wavelength 



modes of the fields are effectively overdamped, and suffer critical slowing down ||3^. Then 
our equations correspond to Model A, in the classification of dynamical stochastic theories 
of Hohenberg and Halperin ||3^. Although Model A will be adequate for our purposes, a 
complete analysis of the universal 0(4) dynamics in the chiral limit of QCD requires coupling 
the order parameter (which is not conserved) to other conserved quantities in the theory, as 
in Hohenberg and Halperin's Model G fl^ . 

For any nonzero value of r], at late times {t ^ f]^^) the Langevin dynamics describes 
evolution towards a steady state characterized by the equilibrium partition function 

Z = N j D(t)E-^'^^*\ (5) 

with 



(6) 



Eq. d^) is the canonical partition function for the linear sigma model in 3D. Note that the 
ensemble of equilibrium configurations described by this partition function is independent of 
rj. In this section and the next, we shall use equilibrium physics to fix all parameters in the 
model except rj. In Section IV, we shall fix 77 by analyzing the dynamics of the relaxation 
of nonequilibrium configurations towards equilibrium. 

Eq. (^ is the correct equilibrium state for the classical model. When compared to the 
quantum partition function it fails to account for the correct statistical weighting of hard 
particles (with typical momentum A; ~ T) and at low temperatures, when (T < m7r,mo-). 
We will return to the latter point in the next section. The classical statistical weighting of 
hard modes in the equilibrium state leads to ultraviolet divergences (as the lattice spacing 
Qs is taken to zero) of certain expectation values at nonzero temperature. This effect is an 
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FIG. 1. The two perturbative diagrams that are ultraviolet divergent in the classical theory at 
finite temperature: a) the tadpole and b) the sunset. Both divergences are momentum independent 
allowing us to compute the sunset diagram at zero external momentum only. 



incorrect prediction of the classical partition function. Fortunately it can be compensated 
for. The idea is to subtract the divergent terms, which are all perturbative and can therefore 
be identified diagrammatically. These are then replaced by the corresponding correct contri- 
butions, identified through the computation of the same diagram using quantum statistical 
distributions. 

For the linear sigma model in three dimensions, there are only two divergent diagrams, 
shown in Fig. |I[ These are self-energy contributions: the tadpole diagram at one loop and 



the sunset diagram at two loops |^5[. Both divergences are momentum independent and can 
be subtracted away by a simple temperature dependent mass renormalization. These two 
diagrams can be written as 



' tadpole 




GM, (7) 



' ^ ' 3 J (2vr)6 |fc + g-12 + ^2 • 

We have made explicit all symmetry factors. The two-point function for the scalar field G{uj) 
has the thermal (no vacuum contributions) form G{uj) = n{uj)/uj, with the single particle 
number distribution given by 



T 

OJ 



classical , 

(9) 

Bose-Einstein , 



for classical or quantum scalar fields respectively. Clearly, the quantum distribution ap- 
proaches the classical for huo/T 0, i.e., in the limit of high occupation numbers. Ulti- 
mately this is why the classical theory becomes a good description of long wavelength physics 
when masses are vanishing or at least small compared to T, as happens at a true critical 
point. 

On the lattice, the continuum dispersion relation for a particle of mass m, namely u"^ = 
k"^ + m^, becomes 

co' = -,j:sm4^]+m\ (10) 
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where is the lattice spacing and D the number of spatial dimensions. Here, ki = 
27mi/Nias, with Ni the linear size of the lattice and rij G { — |A^j, |iVj}, for each Euclidean 
dimension. Eq. (|I0|) reduces to the continuum result when ki is small. 

If we are interested only in the ultraviolet behavior of the tadpole diagram Fig. la, we 
can neglect the mass m. (This diagram has no infrared divergences.) This yields 



/.t„.„^^ATA^0.25(™, (U) 



N + 2XT'^ 

12 h 



The classical tadpole is linearly divergent with the ultraviolet momentum cutoff A. The 
arrow shows the result on the lattice, from summing over modes within a Brillouin zone. To 
the accuracy quoted, these lattice results are identical for Ni > 32 and various temperatures. 
The Bose-Einstein tadpole is finite. 

The sunset diagram Fig. |I]b, has only a logarithmic divergence and so in addition to 
being ultraviolet divergent it is infrared divergent in the absence of masses ^|. Since 



only the difference of the two contributions /^~^ — = AI enters our calculation, the 
infrared divergence cancels out (allowing us again to neglect the mass m) and the ultraviolet 
divergence is (A^ + 2)A^T^ ln(T//iA)/87r^. There is no convenient analytic form for the finite 
part, but on the lattice the result is 

A/sunset ^ [0.0144 ln(ra,/^) - 0.0369] (A^ + 2)\^T^ . (13) 

The renormalization is then achieved by making the replacement 

- XV^ -XV^ + A/tadpole + A/,unset (14) 

in the term linear in (pa in the equation of motion Eq. (0), which removes the linear and 
logarithmic divergences in the self energy. Note that AI as T — >■ 0, and so this renor- 
malization does not affect the choice of the parameters in Eq. (|^) at zero temperature. This 
is true to all loop orders. The procedure we have described implements the matching of 
the classical thermal theory to the quantum in the ultraviolet, where the classical partition 



function explicitly leads to divergences. After adding the counterterms (|Tj); the long wave- 
length classical physics is non-divergent as the lattice spacing is taken to zero. Furthermore, 
these counterterms alone ensure that at high temperature, the long wavelength physics is 
correct and independent of the ultraviolet cutoff. 

At low temperatures, however, differences between the classical and quantum theories 
which do not diverge with the ultraviolet cutoff become important. Such ultraviolet-finite T- 
dependent counterterms arise from many diagrams, not just from the two in Fig. 1. We will 
deal with this pragmatically, by requiring that the behavior of our partition function matches 
the predictions of chiral perturbation theory for T ^ and shows a critical temperature 
Tc in agreement with that found in lattice calculations of full QCD. In the next section, we 
shall use these two criteria to fix the coefficients of finite counterterms which are linear and 
quadratic in T. 

For numerical work, it is convenient to remove the large coupling A ~ 20 from the 
Langevin evolution and express all quantities in dimensionless units by scaling to new primed 
quantities: 
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X 



y/Xvx, t' = ^/\vt, (15) 



H 

V ' Xv 



H'=— (16) 



n = -7^^ T = T. (17) 

Using Eq. (^ gives H' = 0.1263. We must be careful with these rescahngs since they are 
inconsistent with the conventional choice h = 1. The action scales as S" — > S'/X and so in 
order to leave the path integral invariant, Planck's constant h must also be scaled: 

h' = Xh . (18) 

The tadpole and sunset counterterms both scale like 

(A/)'^^. (19) 

Note that we do set /c^ = c = 1 throughout. 

We solve the Langevin equation on a cubic lattice using spatially periodic boundary 
conditions, lattice spacing a^, and time discretization with time step a^. We shall compare 
results obtained using = 0.4 and = 1, in dimensionless units. We always use at = 0.01. 
On the lattice, Eq. (^) is satisfied by choosing ^ independently at each discrete point in space 
and time from a random distribution with (^) = and (^^) = 2ri{T)T / {a^at) . We evolve 
the equations of motion using the fourth order symplectic algorithm, although we have also 
verified that using the second order algorithm suffices. We have found that these algorithms 
introduce much less Oi-dependence than the staggered leapfrog or stochastic Runge-Kutta 
algorithms. 

III. ORDER PARAMETER AND CORRELATION LENGTHS 

To analyze the transition we have measured the order parameter (0) at several temper- 
atures, as shown in Fig. 0. We define (0) by 




(<^) = ^ E h^ E U^.J.k) , (20) 



which is analogous to the magnetization in spin systems. Error bars in the figure are the 
standard deviation from the mean in an ensemble of a few hundred measurements along a 
thermalized Langevin trajectory. To determine after what time the configuration has ther- 
malized, we monitor the system for equipartition and follow the order parameter until well 
after its behavior becomes non-monotonic.[| For H = the system suffers from critical slow- 



To check for equipartition, we only perform the simple test of verifying that the spatial average 
of all four (dcpa/dt)'^ are (like the order parameter) sufficiently time-independent and are all given 
by T. This is equivalent to verifying that the kinetic energy of each field component is T/2. In 
effect, we are measuring the temperature of the fields, which we make sure agrees with the input 
value. 
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FIG. 2. The order parameter (0) as a function of temperature T, for if = (circles) and H 
nonzero (triangles). Open and filled symbols show results obtained with = 1 on a 26^ lattice 
and tts = 0.4 on a 64^ lattice respectively, with counterterms as described in the text. The line 
shows the prediction from chiral perturbation theory, Eq. (21). 



ing down near its critical temperature, and long Langevin trajectories are needed in order 
to obtain thermalization and a statistically meaningful set of measurements. For H 0, 
thermalization occurs quickly at all temperatures, within a time t of order (10 — 100)?7~^. 
Langevin algorithms may yield a^-dependent results if at is not chosen small enough. To ex- 
clude this possibility, we checked our order parameter measurements against those obtained 
from configurations thermalized using a standard Metropolis Monte Carlo algorithm. 

In Fig. H, we have required that the behavior of (0) agrees with the expectations from 
chiral perturbation theory at low temperatures [p7[]. 



(0) = /J 1 



12 f I 



(21) 



In particular, we have enforced the absence of linear T-dependence at small T. As anticipated 
in the previous section, this requires the addition of a further mass counterterm. This 
counterterm is not divergent in the — limit, but it does vary with a^,. It should 
be thought of as coming from the order with p > contributions of infinitely many 
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diagrams. We find that the hnear T-dependence at small T is removed upon making the 
replacement — Af ^ — Af ^ + 6iT, with hi = 0.425 in dimensionless units. Next, we have 
fixed the finite counterterm proportional to T^. In principle, we could have done this by 
enforcing agreement with the coefficient of the term in (|2l|) . This is difficult in practice. 
Instead, we note that in the absence of any finite counterterm proportional to the second 
order phase transition (for H = 0) occurs at Tc ^ 130 MeV, which is somewhat lower than 
that expected in QCD PRp2|. We have pushed Tc up to Tc ~ 150 MeV by introducing 



— Af ^ — >• — Af ^ + biT + 62^^, with 62 = —0.066 in dimensionless units. These values of bi 
and 62 were obtained with a lattice spacing = 0.4 on a 64^ lattice, as shown by the filled 
symbols in Fig. 2. To give the reader a sense of how they depend on a^, note that on a 
26^ lattice with = 1 (and thus the same physical volume as above) we find that with 
62 = 0.084 and hi unchanged from above, the order parameter as a function of temperature 
is the same as that for = 0.4 within error bars, as shown by the open symbols in Fig. 
2. We have verified that the a^-dependence of 62 vanishes in the weak-coupling limit and is 
non-divergent in the ^ limit. 

With H = 0, the order parameter shows a relatively sharp phase transition at T = Tc ~ 
150 MeV, which becomes sharper as the volume of the lattice is increased. At nonzero H, the 
explicit symmetry breaking favors the sigma direction {(pa) = Soafn and the phase transition 
is smoothed into a crossover, which occurs over a range of temperatures somewhat above 
the Tc for H = 0. We shall see below that this crossover is centered at Tcross — 180 MeV. 

In addition to measuring the order parameter we have also measured the temperature 
dependence of the masses of the pion and sigma fields. To achieve this, we introduce a small 
amplitude displacement away from equilibrium in either the sigma or one of the pion field 
directions, and time the resulting oscillations. In practice, after first thermalizing the fields 
at a desired temperature as described above we add an offset in the value of the pion or sigma 
field to the thermalized field configuration, choosing the same offset at every point in space. 
We typically shift one of the pion fields by up to O.l/^r or shift the sigma field by up to 0.05/^. 
(A larger displacement of the sigma leads to non-linear response.) We repeat this procedure 
for ten different realizations of a canonical thermal ensemble. The time-dependent response 
to the perturbation away from equilibrium is most easily measured after first averaging 
(7r)(t) or (cr)(t) over the ensemble of runs. For small enough displacements we expect the 
response to be linear, and characteristic of a damped harmonic oscillator. Thus 

(0a) = (0a(t = 0)) exp[-t/r] cos(u;t), (22) 



with u = \Jm?(T) — t""^ for a spatially homogeneous perturbation. Taking 77 small, in 
order to obtain large relaxation times r, we have measured this behavior at many different 
temperatures. An example of the data and the fit to the form ( [2^ ) is shown in Fig. ^ The 
results we obtain by following this procedure at many temperatures are shown in Fig. ^, 
together with the chiral perturbation theory prediction for the pion mass 



<(r) = ^^(0)(l + ^) , T«/_ (23) 

We thus confirm that with our choice of counterterms, we are doing a good job of reproducing 
the predictions of chiral perturbation theory at low temperatures. We see that the sigma 



10 



1 1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

2 4 6 8 10 12 14 16 18 20 



t[fm] 



FIG. 3. An example of the fits used to determine niT^ and ttt-o- as well as the associated decay 
rates. An instantaneous small perturbation away from the thermal expectation values of the fields 
is induced and its subsequent relaxation is fit to the form (22). In this example, we chose rj = 0.05 
in dimensionless units. 



becomes lighter in the medium as a result of symmetry restoration at high temperatures. 
nicr displays a minimum around Tcross — 180 MeV, which can be taken to be the crossover 
temperature. At Tcross, ""^tt — 220 MeV and m„ ~ 280 MeV. At higher temperatures, mo- 
grows and continues to grow monotonically. The two masses are equal within error bars 
for T > 220 MeV. 

If nature were closer to the chiral limit, m^r and would be less than T for a range 
of temperatures near Tcross, justifying our use of classical field theory for dynamics below. 
As it is, we shall be applying these methods somewhat beyond their regime of quantitative 
validity. 



IV. DISORIENTED CHIRAL CONDENSATE DECAY 



We have now characterized all the parameters and counterterms in the model except 
ri{T). 7] plays no role in the equilibrium physics, but it is crucial to the dynamics. In near- 
equilibrium dynamics, rj is the parameter which controls the rate at which equilibrium is 
approached. To complete the characterization of our linear sigma model, we investigate the 
linear response timescale r for the decay of a long wavelength pion perturbation, excited as 
described in the previous section. For each temperature T, we extract r by fitting the DCC 
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FIG. 4. The thermal masses of the a and vr fields as functions of temperature. Error bars 
represent statistical uncertainty in the ensemble of measurements and uncertainty in the fit, see 
Fig]^. The line shows the prediction from chiral perturbation theory, Eq. ([2^). The minimum of 



rUn 



which is a reasonable definition of the crossover temperature, is reached around 



180 



MeV. All calculations were done with a, = 0.4 on a 64^ lattice. 



decay profile to the form (p2D . The r so obtained depends on our choice of input i].^ 

The timescale t(T) for the decay of a region of disoriented chiral condensate in contact 



with an equilibrated gas of hadrons has been computed previously by Steele and Koch [19 



We shall seek to choose ri{T) so that the t(T) we then extract reproduces that calculated 
as in Ref. 



19 



In so doing, we are normalizing the magnitude of the dissipation in the 
Langevin evolution to a physical observable. 

The dissipation parameter r]{T) has been computed perturbatively by Rischke 



20| from 

the imaginary part of the two-loop pion self-energy in the finite temperature quantum linear 
sigma model, with the result: 



^Note that we only measure r for the decay of a pion perturbation. It may also be interesting 
to investigate the timescale for the relaxation of a perturbation in the sigma direction, r and 
To- will differ at low temperatures and become degenerate above Tcross- In such an investigation, 
one could introduce rja-a 7^ rjirn, and study how r and Tq- depend on r/o-o- and ry^yr- Note that even 
though in this paper we take r^o-o- = Vim = Vi this does not mean that r and are equal. 
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2 X — e-^^/^f^TT g(m2-2m2)/2m^T _ ^ ' 



(24) 



At A = 20, we can only expect this calculation to be a qualitative guide. 

We shall slightly extend the calculation of t(T) performed by Steele and Koch, and so 
instead of simply quoting the result we first sketch their computation. We wish to analyze a 
large, spatially uniform, DCC in contact with a thermal gas of many other hadronic species 
X = n, K, which will interact with the DCC and ultimately cause it to decay. Modeling 
the DCC by a coherent state, and working through the LSZ formalism, Steele and Koch 
show that the depletion of the particles in the DCC, dN/dt, is proportional to the square of 
the interaction amphtudes |19|: 



dN 
It 



d^ki d^k2 d^k^ 



(27r)3 (27r)3 (27r) 



IC 



X 



fc2+fc3 — fcl 

i2Eoy 



(25) 



For TT-vr scattering, the thermal weighting is given by 



123 



/2/3(l+/l)-/l(l + /2)(l + /3 



(26) 



1^ 



^ representing the Bose-Einstein momentum distribution. The 
[{k2 + k^ — kiY + m^]-^/^. \C\'^ to first approximation gives the 



with fi = {exp{Ei/T) - 
DCC pion energy is Eq 
momentum-conserving delta-function leading to a formula for the half-life r of the DCC 



1 dN 



d^ki d^k2 d^ks 



N dt 2EoJ (27r)3 (27r)3 (27r) 



X 



X 
123 



r 



(27) 



The focus of Ref. [|rU| was on the late time decay of a DCC, due to a heat bath with a 
temperature well below Tc. In this context it is reasonable to take the scattering amplitudes 
T^x from phase-shifts of real zero-temperature data. If these are not known, then a Breit- 
Wigner form can be fit to the dominant resonance. To obtain the values of t(T) shown below, 
TTTT was taken from data, tiK was modeled by the ii'*(892) resonance, vrp was modeled by 
the ai resonance, ttN was taken from data, TriV was taken from the same vrA^ data, and 
ttA (ttA) was modeled by the iV*(1675) resonance. Note that at T = 170 MeV the pions 
account for only about half the particles in our equilibrated resonance gas. All the other 
species together do play an important role. We have extended the previous analysis by 
including the dissipation due to ttK, ttA and vriV scattering, not included in Ref. [|T^]. The 
r(T) obtained in this way is shown in Fig. ^ The calculation is only quantitatively valid 
below Tc, but we have plotted the results up to higher temperatures also. 

Below we shall choose ti{T) to match the decay time t(T) of a spatially homogeneous 
DCC in our Langevin evolution to the values calculated by Steele and Koch. First, however, 
let us see what we obtain with rj = 0, as in Fig. ^. That is, we first thermalize at a 
given temperature, add a DCC, and then set r] = thus removing the classical fields from 
further contact with any heat bath. In so doing, however, we do not remove all sources 
of dissipation. The long wavelength DCC still decays, by virtue of its couphng to those 
finite wavelength modes which we are treating classically. That is, the DCC decays via 
contact with an effective classical heat bath, present even for rj = 0. We see that at low 
temperatures the resulting Tci(T) turns out to be smaller than the desired t(T), namely that 
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20 40 60 80 100 120 140 160 180 200 220 

T [MeV] 



FIG. 5. The solid curve shows the decay time of a spatially homogeneous DCC via interaction 
with a gas of hadrons, computed following Ref. |jl^ as described in the text. The open (filled) 
circles show Tci(T), obtained from the Langevin equation with rj = 0, for as = 1 {as = 0.4). 
The uncertainty in each point (statistical and that from the fit) is comparable to the dispersion 
between open and closed circles. Tci(T) reflects the dissipation of the long wavelength DCC induced 
by contact with thermalized classical fields at shorter wavelengths. We see that Tci(T) < t(T) below 
T ~ 145 MeV. 

due to contact with a gas of hadrons. The reason is that t(T) becomes very large at low 
temperatures, due to the exponential suppression of the thermal occupation numbers in the 
(quantum) gas of hadrons. The classical thermal field theory overestimates the occupation 
numbers, and thus the dissipation, leading to an underestimate of r. 

Note that we have checked that our calculation of Td is independent of the lattice spacing 
ag. As shown in Fig. 5, we have compared results obtained on a 64^ lattice with = 0.4 
to those obtained on a 26^ lattice with = 1. As described in the previous section, 
the ultraviolet-finite counterterms must be chosen differently in these two cases in order 
that the equilibrium behavior of the order parameter is the same in both cases. With the 
counterterms so adjusted, TdiT) changes little with a^. 

The shape of Tci(T), shown in Fig. ^, is interesting and may be understood as follows. 
Tci{T) should grow like 1/T as T ^ (whereas r grows much more rapidly in this limit). 
We do find that Tc\ rises at the lowest temperatures we have explored. The minimum of Td 
coincides with the minimum value of the ratio m.„^{T)/T, which occurs at T ~ 120 — 130 MeV. 
This means that the classical thermal pions are most numerous in this temperature regime, 
making it plausible that the classical dissipation is largest and Td is smallest. 
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FIG. 6. Circles show the DCC decay time Tci(T) for the classical u-model at r] = 0. Squares 
show t(T) obtained from the Langevin equation with 7?(T) taken from Rischke's perturbative 
calculation. Diamonds show t{T) obtained from the Langevin equation with ri{T) chosen so as to 
reproduce t(T) from the calculation of Steele and Koch, shown as the solid line. Errors are similar 
to those discussed for Fig. 5. 



We now see the challenge posed by our use of a classical theory, with its intrinsic decay 
time Tci{T). By turning on a nonzero f], we can only add dissipation. This means that we 
can only reduce r, relative to Tci{T). At temperatures where Tci(T) < r(T), therefore, there 
is no way to use our Langevin equation to reproduce the correct dynamics of the decay of 
a DCC. Thus, we find that our model can only be used to describe the long wavelength 
dynamics in the presence of a heat bath with T ^ 145 MeV, corresponding to about 80% 
of the crossover temperature Tcross- T ~ 145 MeV is large enough that the assumptions in 
the calculation of r as in Ref. [T^ may be starting to break down, which means that our 
estimate of the limit of validity of our analysis may have a little play in it. However, the 
fundamental fact that r should rise rapidly with decreasing temperature while Td does not 
is incontrovertible. 

Above T ~ 145 MeV, we can proceed as planned. We can choose an f]{T) in such a way 
as to reproduce the predictions of Steele and Koch, as demonstrated in Fig. |. The required 
ri(T) is plotted in Fig. 0, where it is also compared to the ri{T) of Eq. (p4D, calculated 
perturbatively by Rischke pO|. We see (from both Figs. ^ and ^ that introducing the 
perturbatively calculated ?7(T) into the Langevin equation is not sufficient to obtain a DCC 
decay timescale t(T) in agreement with that induced via contact with a thermal hadron 
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FIG. 7. ri{T) from Rischke's perturbative calculation (solid line) and r](T) chosen so that the 
Langevin equation reproduces t(T) computed by Steele and Koch (circles). 



gas.| Both the calculations of Rischke and of Steele and Koch are only of quantitative 
validity below T^, as they use zero temperature masses and interactions. Nevertheless, the 
discrepancy between the ?7(T)'s in the vicinity of is sufficiently large that we take this to 
mean that, for A = 20, the perturbative calculation underpredicts rj. 



V. DISCUSSION 

We have shown how to use a stochastic classical linear sigma model to obtain a reasonable 
description of the equilibrium thermodynamics and the long-wavelength near-equilibrium 
dynamics of QCD at temperatures of order that of its phase transition. The use of any 
classical statistical field theory around the chiral transition cannot be quantitatively jus- 
tified, because there are no excitations with masses which are very small compared to the 
temperature. Nevertheless, we have determined the choice of divergent counterterms needed 



^It is not clear whether we should have expected the use of t] from ( p^ ) to yield too rapid or 
too slow dissipation. On the one hand, we may have expected that use of the perturbative 7]{T) 
in the Langevin equation would yield excess dissipation, as there is always classical dissipation in 
addition. On the other hand, we may have expected that it would introduce insufficient dissipation, 
as its calculation neglects the effects of all hadrons in the heat bath except the pions and sigma. 



16 



to remove the ultraviolet divergences from the classical linear sigma model, making the long 
distance physics described therein independent of the lattice spacing. We have then been 
able to find the finite counterterms needed in order for the equilibrium temperature depen- 
dence of the order parameter in the linear sigma model to be as in QCD at low temperatures, 
as calculated in chiral perturbation theory, and for Tc to be as in QCD, as calculated in lat- 
tice gauge theory. The temperature dependence of the pion and sigma masses then come 
out nicely in accord with expectations. We have then fixed rj{T) such that the timescale 
for the dissipation of a DCC in our Langevin model agrees with that calculated previously. 
The main constraint on our model is that it cannot be used at low temperatures, where 
the classical dissipation is too great even when rj is set to zero. We find that this limits 
its use to temperatures above about 80% that of the crossover transition itself. Above this 
temperature, we show how to choose ?7(T) in such a way that the long wavelength dynamics 
of the pions and sigma should be reproduced reasonably well in our model, at least as long 
as the system is not too far out of thermal equilibrium. 

In this work, we have used the dynamics of an idealized DCC as a device with which 
to normalize our model. In subsequent work, we shall use the now normalized model to 
investigate the dynamical setting where the system expands and cools through its chiral 
crossover transition. We hope to determine with some confidence how fast the expansion 
and cooling must be in order for the dynamics to be driven significantly from equilibrium. 
We will then follow the nonequilibrium long wavelength evolution for as long as occupation 
numbers remain reasonably high. Unlike in the idealized setting utilized for convenience 
in this paper and in the analytic computations of rjiT), the long wavelength modes will 
certainly not describe a spatially uniform condensate. 



ACKNOWLEDGEMENTS 

KR is grateful to D. Rischke for long discussions several years ago which started DR down 



the path toward Ref. [20] and started KR down the path toward the present paper. We also 
acknowledge conversations with A. K. Chaudhuri. Numerical work was performed in part 
at the T-Division/CNLS Avalon cluster at Los Alamos National Laboratory. This work 
is supported in part by the Department of Energy under cooperative research agreement 
#DF-FC02-94ER40818. The work of KR was supported in part by a DOE OJI Award and 
by the A. P. Sloan Foundation. 



17 



REFERENCES 



[1] R. Baier, A. H. Mueller, D. Schiff, and D. T. Son, Phys. Lett. B 502, 51 (2001) jh^ 
ph/0009237| . 



[2] K. H. Ackermann et al. [STAR Collaboration], Phys. Rev. Lett. 86, 402 (2001) inucT 
ex/00090ll| 



[3] A. Barducci, R. Casalbuoni, S. DeCurtis, R. Gatto, G. Pettini, Phys. Lett. B 231, 463 
(1989); S.P. Klevansky, Rev. Mod. Phys. 64, 649 (1992); A. Barducci, R. Casalbuoni, 
G. Pettini and R. Gatto, Phys. Rev. D 49, 426 (1994); M. Stephanov, Phys. Rev. 
Lett. 76, 4472 (1996); Nucl. Phys. Proc. Suppl. 53, 469 (1997); J. Berges and K. Ra- 
jagopal, Nucl. Phys. B 538, 215 (1999) ||hep-ph/9804"233i ; M. A. Halasz, A. D. Jackson, 
R. E. Shrock, M. A. Stephanov and J. J. Verbaarschot, Phys. Rev. D 58, 096007 (1998) 
||hep-ph/ 98042901 . 

M. Stephanov, K. Rajagopal and E. Shuryak, Phys. Rev. Lett. 81, 4816 (1998) ^epj 
ph/98062l9| . 



[6] 
[7] 



[9 
[10 

[11 

[12 
[13 



[14 

[15 



M. Stephanov, K. Rajagopal and E. Shuryak, Phys. Rev. D 60, 114028 (1999) |iep^ 
ph/9903292|| . 



For a review see K. Rajagopal, Acta Phys. Polon. B 31, 3021 (2000) [|hep-ph/0009058 
K. Rajagopal, Nucl. Phys. A 680, 211 (2000) [ |hep-ph/0005T0T| . 

See, for example, H. Heiselberg and A. D. Jackson, |nucl-th/9809013| ; I. N. Mishustin, 
Phys. Rev. Lett. 82, 4779 (1999) ||hep-ph/9811307| . 

B. Berdnikov and K. Rajagopal, Phys. Rev. D 61, 105017 (2000) ||hep-ph/9912"27i . 



K. Rajagopal and F. Wilczek, Nucl. Phys. B 404, 577 (1993) ||hep-ph/ 93032811 . 
A. A. Anselm, Phys. Lett. B 217, 169 (1988); A. A. Anselm and M. G. Ryskin, Phys. 
Lett. B 266, 482 (1991). 

J. Blaizot and A. Krzywicki, Phys. Rev. D 46, 246 (1992). 

J. D. Bjorken, Int. J. Mod. Phys. A7, 246 (1992); ibid., 4189; J. D. Bjorken, Acta 
Phys. Polon. B 23, 637 (1992); K. L. Kowalski and C. C. Taylor, |hep-ph/921128^ ; 
J. D. Bjorken, K. L. Kowalski and C. C. Taylor, |hep-ph/9309235| . 
K. Rajagopal and F. Wilczek, Nucl. Phys. B 399, 395 (1993) ||hep-ph/921025l . 
M. M. Aggarwal et al. [WA98 Collaboration], Phys. Lett. B 420, 169 (1998) |Eepj 
I ex/9710"0T5| . M. M. Aggarwal et al. [WA98 Collaboration], Phys. Rev. C 64, 011901 

(2000) ||nucl-ex/00l2001 . 

[16] C. Greiner, C. Gong and B. Muller, Phys. Lett. B 316, 226 (1993) [|hep-ph/9307336 



Z. Huang, M. Suzuki and X. Wang, Phys. Rev. D 50, 2277 (1994) | |hep-ph/9403'5UD 
S. Gavin, Nucl. Phys. A 590, 163 (1995); |hep-ph/9407368| ; Z. Huang, I. Sarcevic, 
R. Thews and X. Wang, Phys. Rev. D 54, 750 (1996) ||hep-ph/9511387|| ; Z. Huang and 
X. Wang, Phys. Lett. B 383, 457 (1996) ||hep-ph/9604300| ; D. Boyanovsky, H. J. de 
Vega, R. Holman and S. Prem Kumar, Phys. Rev. D 56, 5233 (1997) ||hep-ph/970T360 



Phys. Rev. D 56, 3929 (1997) ||hep-ph/ 97034221 ; Y. Kluger, V. Koch, J. Randrup and 
X. Wang, Phys. Rev. C 57, 280 (1998) [|nucl-th/9704018|| ; H. Hiro-Oka and H. Minakata, 
Phys. Lett. B 425, 129 (1998) [Erratum-ibid. B 434, 461 (1998)] [|hep-ph/9712476 



T. C. Petersen and J. Randrup, Phys. Rev. C 61, 024906 (2000) [|nucl-th/9907051 



C. Chow and T. D. Cohen, Phys. Rev. C 60, 054902 (1999) | |nucl-th/9908TJT3| ; J. I. Ka- 
pusta and S. M. Wong, Phys. Rev. Lett. 86, 4251 (2001) ||nucl-th/001200^] . 
[17] F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz and P. R. Anderson, Phys. Rev. 



18 



D 50, 2848 (1994) ||hep-ph/9405352|; D. Boya novsky, H. J. de Vega, and R. Holman, 
Phys. Rev. D 51, 734 (1995) ||hep-ph/ 940 13081 ; F. Cooper, Y. Kluger, E. Mottola, and 
J. Paz, Phys. Rev. D 51, 2377 (1995) ||hep-pli/9401308| ; M. A. Lampert, J. F. Dawson 
and F. Cooper, Phys. Rev. D 54, 2213 (1996) ||hep-th/960306g] ; F. Cooper, Y. Kluger, 
E. Mottola, Phys. Rev. C 54, 3298 (1996) ||hep-ph/ 9604281 ; D. Boyanovsky, F. Cooper, 
H. J. de Vega and P. Sodano, Phys. Rev. D 58, 025007 (1998) ||hep-ph/ 980227^ . 
[18] M. Asakawa, Z. Huang and X. Wang, Phys. Rev. Lett. 74, 3126 (1995) | |hep-ph/940g2gU 



[19] J.V. Steele, and V. Koch, Phys. Rev. Lett. 81, 4096 (1998) [^lucl-th/ 9806055 1. 
[20] D. H. Rischke, Phys. Rev. C 58, 2331 (1998) 



nucl-th/ 98060451]. 



[21] C. Greiner and B. Muller, Phys. Rev. D 55, 1026 (1997) [ |hep-th/ 9605048 



[22] Z. Xu and C. Greiner, Phys. Rev. D 62 036012 (2000) [|hep-ph/9910562 



[23] T. S. Biro and C. Greiner, Phys. Rev. Lett. 79, 3138 (1997) ||tiep-ph/ 97042501 . 
[24] A. K. Chaudhuri, Phys. Rev. D 59, 117503 (1999) ||nucl-th/98090T^ ; J. Phys. G 27 175 
(2001) ||hep-ph/990837'B| ; ^ep-ph/0007331 ; |hep-ph/OOllOOl . 



[25] A. K. Chaudhuri, |hep-ph/01050T3 



[26] J. Randrup, Phys. Rev. Lett. 77, 1226 (1996) [|hep-ph/9605223|| ; Nucl. Phys. A 616, 
531 (1997) ||hep-ph/ 96 12451 ; J- Randrup, Phys. Rev. C 62, 064905 (2000) fmcT- 

I th/ooioioi . 

[27] R. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984); F. Wilczek, Int. J. Mod. Phys. 
A 7, 3911 (1992). 

[28] S. Gottlieb et al., Phys. Rev. D 55, 6852 (1997) ||hep-lat/ 96 12020 



[29] R. D. Pisarski, Phys. Rev. D 62, 111501 (2000) ||hep-ph/0006"205|] ; A. Dumitru and 

R. D. Pisarski, Phys. Lett. B 504, 282 (2001) ||hep-ph/0010"083| . 
[30] S. Datta and S. Gupta, Nucl. Phys. B 534, 392 (1998) ||hep-lat/ 9806031 . 
[31] R. V. Gavai and S. Gupta, Phys. Rev. Lett. 85, 2068 (2000) ||hep-lat/0004UTT 



[32] For reviews, see F. Karsch, Nucl. Phys. Proc. Suppl. 83, 14 (2000) [hep-lat/9909006|| ; E. 

Laermann, Nucl. Phys. Proc. Suppl. 63, 114 (1998); and A. Ukawa, Nucl. Phys. Proc. 

Suppl. 53, 106 (1997). 

[33] L. M. A. Bettencourt, Phys. Rev. D 63, 045020 (2001) ||hep-ph/000526l . 

[34] P.C. Hohenberg and B.L Halperin, Rev. Mod. Phys. 49, 435 (1977). 

[35] G. Parisi, Statistical Field Theory (Addison- Wesley, New York, 1988). 

[36] K. Farakos, K. Kajantie, K. Rummukainen, M. Shaposhnikov, Nucl. Phys. B 425, 67 

(1994) ||hep-ph/940420T| . 
[37] P. Gerber and H. Leutwyler, Nucl. Phys. B 321, 387 (1989); R. D. Pisarski and M. Tyt- 

gat, Phys. Rev. D 54, 2989 (1996) ||hep-ph/ 9604401 . 



19 



